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On  the  uniformly  minimum  variance  unbiased  estimators 


of  the  variance  and  its  reciprocal  of 
an  Inverse  Gaussian  distribution 

Abstract 

The  two-parameter  Inverse  Gaussian  distribution  is  found  to  have 
useful  applications  in  a wide  variety  of  fields.  The  uniformly 
minimum  variance  estimator  of  its  mean  is  known  and  is  the  sample 
mean.  However,  no  such  estimator  of  the  variance  is  reported  in 
literature.  Here  the  uniformly  minimum  variance  unbiased  estimators 
of  the  variance  and  its  reciprocal  are  derived. 


KEY  WORDS:  Inverse  Gaussian  distribution;  reciprocal  of  the  variance; 
variance;  uniformly  minimum  variance  unbiased  estimator. 
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On  the  uniformly  minimum  variance  unbiased  estimators 
of  the  variance  and  its  reciprocal  of 
an  Inverse  Gaussian  distribution 

1.  INTRODUCTION 

The  Inverse  Gaussian  distribution  was  introduced  by  Schrbdinger 
(1915),  and  has  been  studied  extensively  by  many  authors, pioneered 
by  Tweedie  (1945,  1956,  1957a,  1957b),  and  Wald  (1945).  The  distri- 
bution is  found  to  have  useful  applications  in  a wide  variety  of 
fields,  such  as  biology,  economics,  medicine,  and  reliability  and  life 
testing.  See  the  recent  review  article  by  Folks  and  Chhikara  (1978) 
and  discussion  on  it. 

Using  Tweedie's  notation  the  probability  density  function  of  an 
Inverse  Gaussian  random  variable  X is  given  by 


f(x;  p,A)  = (A/2irx^)2  exp  (-A(x-p) 2/2p^  x},  x>0, 

= 0 , otherwise 


(1.1) 


where  p,A  are  both  positive.  As  far  as  estimation  of  the  parameters 
p and  A is  concerned,  it  is  known  that  the  sample  mean  is  the 
uniformly  minimum  variance  unbiased  estimator  of  the  mean  p of  X. 
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However,  no  such  estimator  has  been  obtained  in  literature  for 
3 

the  variance  y /X  . See  Folks  and  Chhikara  (1978).  Theorem  1 

below  gives  this  estimator,  and  Theorem  2 gives  the  uniformly 

. • 3 

minimum  variance  unbiased  estimator  of  the  reciprocal  A/y  of 

the  variance. 

2.  THE  ESTIMATORS  AND  THEIR  DERIVATIONS 
Let  X|,...,Xn(n  > 2)  be  a random  sample  form  ( 1 . 1 ) , and  X 

n i i 1 

be  the  sample  mean,  and  v = E Y ~ = • Then 

i=l  *i  x , 

Theorem  1.  The  uniformly  minimum  variance  unbiased  estimator  of 
the  variance  of  (1.1)  is  given  by  for  n > 3 

T(Xr...Xn)  = nX2  (1  + U)  (2.1) 

where  U is  given  by 


where  k = (n-3)/2  if  n is  odd  = (n-4)/2  if  n even,  and 

W = 2tan-1  /(VX)/n  if  n is  even  and  = n if  n is  odd,  and 
where  B( p , q)  is  the  usual  beta  integral.  And  for  n = 2 the 
estimator  is  given 


T(x  x ) = 

1 ‘ VX  + 2 


Pxoo fj:  It  is  known  that  the  statistic  (X,  V)  is  complete  and 

o n P 

sufficient  for  the  family  (1.1).  The  sample  variance  S = 7.  (X.-X)  /( n-1 ) 

1 ' 

2 3 

is  unbiased  for  the  variance  a = y /k  . To  finish  the  proof  for 

o 

the  case  n = 2,  we  apply  the  Lehmann-Scheffe  Theorem  to  S and 
note  that  S2  = 2 X 3V/(VX  + 2)  . 

The  proof  for  the  case  n >_  3 is  however  more  involved. 

For  an  application  of  the  Lehmann-Scheffe  Theorem  using  $/_,  we  need  to 

2 

compute  the  conditional  expectation  of  X^  given  (X,  V)  . The 
conditional  density  hlx^lc.v)  of  X-j  given  (X,  V)  is  derived 
in  Chhikara  and  Folks  (1974)  in  a different  context.  It  is  given 
by,  for  n 3 , 


h(x-j  ( x,v) 


/n  (n-1) 

1 n-2 

V 2 

73 


vx3  (nx  - x1 ) 3 


U 


1 - 


J L 


— 2 

n(x1  - x) 

vxx;  ( nx-x, ) 

1 1 

(2.2) 


n-4 


where  L < x-j  < U and  L and  U are  the  two  roots  of  the 

o 

quadratic  equation  n(x-j  - x)  = vx^x(nx  - x-j)  . Now  letting 


u = /n(x-|-x)/{vxx^  (nx  - x-j)}^  in  (2.2)  we  obtain  after 

simpl i fi cations , 


some 


E(X^|x,v)=  nx 


- 2 


2( n-1  )nx  c T 

'T^r  1 

2’  2 


B 


(2.3) 


where 


r 1 


1 


0 n+vxu 


ZTT  (1-u  ) 


n-4 
2\  2 


du 


(2.4) 


V 
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To  evaluate  the  integral  I in  (2.4)  , first  write  the  integrand  as 


_L 

vx 


k 

E 

r=l 


r-1  k_r 

vx  + n 

vx 

r\  K I 

(1-u2) 

Ik 

vx  + n i 


vx 


I 

~2  ,n  even 

n+vxu4 


(2.5) 


_L 

vX 


k 

E 

r-1 


vx  + n 


vx 


r-1 


(1-u2) 


>k 


k r" 2 + (vx  + n } 

{ vx  J 


n-uV 

n+vx  u2 


n odd 


where  k = (n-4)/2  if  n is  even  and  = (n-3)/2  if  n is  odd. 

(In  deriving (2. 3)use  is  made  of  the  fact  that  x^  = nx(2  +vxu  )/2(n+vxu  ) 
+g(u),  where  g is  an  odd  function  of  u.j  Now  to  finish  the  proof 

r, 

use(  2.5)  in  (2.4)and  standard  integrals,  and  note  that  E(S  ;x,v) 

= n E(X2 |x,v)/(n-l)  - n2x2/(n-l)  . 


Theorem  2:  The  uniformly  minimum  variance  unbiased  estimator 
of  the  reciprocal  of  the  variance  of  (1)  is  given  by,  for  n .>  4 , 


T(X] 


v ) = ("-3)  _ JL_  + JL 

n VX J nX  n x 


V 

TT-T)  • 


(2.6) 


Proof:  The  unbiasedness  of (2.6)follows  from  the  expressions 
for  negative  moments  about  zero  of  an  Inverse  Gaussian  distribution 
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/ 


given  by  Tweedie  (1957a)  and  the  facts  that  (i)  J and  V are 
independent  and  that  (ii)  X has  an  Inverse  Gaussian  distribution 
with  parameters  (p,  nA),  (iii)  A V has  a chi-square  distribution 
with  (n-1)  degrees  of  freedom.  An  application  of  the  Lehmann- 
Scheffe  Theorem  to(2.6)using  the  fact  that  (X,  V)  is  complete 
as  well  as  sufficient  finishes  the  proof. 
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